A model for a non-minimally coupled scalar field interacting with dark matter 
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In this work we investigate the evolution of a Universe consisted of a scalar field, a dark matter 
field and non-interacting baryonic matter and radiation. The scalar field, which plays the role of 
dark energy, is non-minimally coupled to space-time curvature, and drives the Universe to a present 
accelerated expansion. The non-relativistic dark matter field interacts directly with the dark energy 
and has a pressure which follows from a thermodynamic theory. We show that this model can 
reproduce the expected behavior of the density parameters, deceleration parameter and luminosity 
distance. 



I. INTRODUCTION 

Recently, well known astronomical observations with 
super-novae of type la suggested that our Universe is 
presently submitted to an accelerated expansion 0, Q 
but the nature of the responsible entity, the so-called dark 
energy, still remains unknown. The simplest theoretical 
explanation for the acceleration consists in introducing 
a cosmological constant and investigating the so-called 
AC DM model (see Q), which fits the present data very 
well but has some important unsolved problems. Another 
possibility is to introduce a minimally coupled scalar field 
4> (t) , which has been extensively studied by the scientific 
community and could be invalidated if future observa- 
tions imply that the ratio between the pressure and the 
energy density of the scalar field is restricted to values 
Pcjy/p^ < -1 (see 01). 

In this work we consider a scalar field non-minimally 
coupled to space-time curvature, which was investigated 
in [5j and widelystudied recently, among others, in the 

works HHHUESmill- Also J-2 side several models for 
the dark sector interaction 0, 0, , we consider here 
a direct coupling between dark matter and dark energy, 
following a model proposed in 0] and studied more re- 
cently in |l7j and 01 • We also use a thermodynamic the- 
ory |l9| in order to relate the effects of the non-minimally 
coupling to the dark matter pressure [l^. All compo- 
nents will be described by a set of field equations, and the 
resulting observables - namely, the density parameters, 
the decelerating parameter and the luminosity distance - 
which arc obtained as solutions of the field equations, will 
be compared to the available data set in order to draw 
the conclusions about the viability of this model. We 
show that physically acceptable solutions are obtained 
and that there exist some freedom parameters that will 
be important to fit the data from incoming experiments. 
Units have been chosen so that 8irG = c = h = k = 1, 
whereas the metric tensor has signature (+,—,—,—). 



II. FIELD EQUATIONS 



and radiation. Furthermore, we shall consider that the 
scalar field, which plays the role of dark energy, is non- 
minimally coupled to the scalar curvature [5j and that 
the dark matter interacts with the scalar field according 
to a model proposed by Wetterich [TjJ C3 • 

For an isotropic, homogeneous and spatially flat Uni- 
verse described by the Robertson- Walker metric 



ds 2 = dt 2 - a{tf (dx 2 + dy 2 + dz 2 ) 



(1) 



with a(t) denoting the cosmic scale factor, the Friedmann 
and acceleration equations read 
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respectively. The quantities p = Pb + p r + Pdm + P4> 
and p — pb + p r + Pdm + P<p denote the energy density 
and the pressure of the sources of the gravitational field, 
respectively, and they are given in terms of the respective 
quantities for its constituents. Moreover, the dot denotes 
the differentiation with respect to time, H — a /a is the 
Hubble parameter, and £ the coupling constant between 
the scalar field <f> and the curvature scalar. 

The modified Klein-Gordon equation for the scalar 
field, which takes into account the non-minimal coupling 
and the interaction with the dark matter, reads 



3H(f> 



dV 
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(p - 3p) - Ppdm = 0, (3) 



where (3 is a constant that couples the fields of dark en- 
ergy and dark matter. The energy density and the pres- 
sure of the dark energy are given by 
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respectively, where V(<fi) denotes the dark energy poten- 
tial density. From - © it follows the evolution equa- 
tion for the energy density of dark energy, namely 



We shall model the Universe as a mixture of a scalar 
field, a dark matter field and non-interacting baryons 



3H (p^ + p^) = - + 0Pdm<]>- (6) 



2 



The first term on the right-hand side of © is the respon- 
sible for the energy transfer from the dark energy to the 
gravitational field whereas the second refers to an energy 
transfer to the dark matter field. 

We shall assume that the baryons are non-relativistic 
particles so that pb = and that the barotropic equation 
of state for the radiation field p r = p r /3 holds. Hence, 
the evolution equations for the energy densities of non- 
interacting baryons and radiation become 



p b + 3Hp b = 0, p r + AHp r = 0, 



(7) 



respectively, and we get from Q the well-known rela- 
tionships: pb oc 1/a 3 and p r oc 1/a 4 . 

Now by using equations 10) together with © and (JTJ 
it follows the evolution equation for the energy density 
of the dark matter field: 



Pdr, 



3H (p d „ 
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(8) 



Note that the term on the right-hand side of (JSJ rep- 
resents the energy transfer from the dark energy to the 
dark matter. 



III. APPLICATION 

Two different situations for the dark matter could be 
analyzed, namely, the one where it is considered as pres- 
sureless and the other with a non- vanishing pressure. 
Here we shall restrict ourselves to the latter case. In 
order to determine the expression for the dark matter 
pressure, we recall that in the presence of matter creation 
the first law of thermodynamics for adiabatic (dQ = 0) 
open systems reads 



dQ = d(pV) +pdV 



Pdr, 



' Pdr, 



Udr 



■d(n dm V) = 0, (9) 



where V denotes the volume and it was supposed that 
only dark matter creation is allowed. By considering 
n-dm oc pdm and V oc a 3 it follows from together 
with © and © - ©: 
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(10) 



Now we use again (JHJ and obtain the following expression 
for the pressure of the dark matter field 



Pdr, 
Pdn 
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iHpdr, 



(ii) 

We note that only the plus sign must be chosen in the 
above equation in order to obtain a positive dark matter 
pressure. Moreover, we observe that for this case the 
dark matter pressure vanishes when f3 = £ = 0. 

We have solved numerically the system of differential 
equations © and JHJ) as functions of the red-shift z = 
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FIG. 1: Evolution of density parameters vs. red-shift. 
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FIG. 2: Deceleration parameter vs. red-shift. 



1/a - 1 by using (JU, ©, (JTT1) together with p b oc 1/a 3 
and p r oc 1/a 4 and the usual exponential potential V = 
Voe~°"^, where a and Vo are constants. Furthermore, 
the initial conditions for the energy densities were chosen 
from the present known values given in the literature 
|20j for the density parameters fli(z) = pi(z)/p(z), i.e., 
h b (0) = 0.04995, 0.(0) = 5 x 10~ 5 and n dm (0) = 0.23. 
The constant Vq was determined from the present value 
for the dark energy density parameter, i.e., 0^(0) = 0.72. 

In Fig. 1 we have plotted the red-shift evolution of the 
density parameters where it was considered a fixed value 
for the exponent a, namely a = 0.1, and different values 
for the coupling constants (3 and £. From this figure 
we infer that all kinds of coupling delay the growth of 
dark matter energy density and the decay of dark energy 
density with respect to the red-shift, by a comparison 
with the uncoupled case, i.e., = ^ = 0. Also, positive 
values of £ imply that the radiation-matter (dark matter 
plus baryons) equality takes place at higher red-shifts, 
whereas negative values of £ at lower ones. 

The deceleration parameter q = —da/ a 2 as a function 
of the red-shift is plotted in Fig. 2. According to cur- 
rent experimental values |2j the decelerated-accelerated 
transition occurs at zt = 0.46±0.13, whereas the present 
value of the deceleration parameter for the AC DM case 
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FIG. 3: Difference vs. red-shift. 
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FIG. 4: Dark matter pressure vs. red-shift. 



is qo « —0.55. We note from this figure that the 
decelerated-accelerated transition shifts from zt ~ 0.7 
for the uncoupled case (/3 = £ = 0) to Zt w 0.2 for 
(3 = —0.5 and £ = 0. Moreover, by keeping a fixed 
value for [3, positive values for £ cause later transitions 
where negative values for £ imply earlier transitions. This 



means that, by adjusting the coupling constants in this 
model, it is possible to have concordance with the exper- 
imental data. 

In Fig. 3 we have plotted the difference no between the 
apparent magnitude m and the absolute magnitude M of 
a source as a function of the red-shift where the circles 
represent the experimental values taken from the work 
by Riess et al. [2l| for 185 data points of super-novae 
of type la. The expression for reads /Lto = m — M = 
5 log cLl + 25, where d^ is the luminosity distance given in 
Mpc: d L = (1 + z)^ 1 J* dz/H(z). We infer from the 
figure that all cases fit well the experimental data at low 
red-shifts, but when £ becomes more negative, the curve 
goes away from the expected experimental values. 

The pressure of the dark matter as function of the red- 
shift is plotted in Fig. 4. We infer from this figure that for 
a fixed value of the coupling constant j3, negative values 
for £ lead to higher pressure of the dark matter field, and 
by increasing £ to positive values the pressure decreases. 
In all cases the dark matter pressure tend to a constant 
value for higher red-shifts. 



IV. CONCLUSIONS 

The model analyzed here includes two kinds of cou- 
plings, namely, the coupling of the scalar field with the 
curvature scalar and the coupling of the dark matter 
field with the scalar field. Like the model investigated in 
the work [TJ, it has a strong dependence on the initial 
conditions, since we have solved the system of coupled 
differential equations by imposing that the slope of the 
scalar field has small positive values at very high red- 
shifts, which was a necessary condition to have positive 
values for the dark energy density. However, as we have 
shown, there exist some freedom in this model because 
the two coupling constants can be adjusted in order to 
fit present and future experimental data. 
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